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RELAXATION  PHENOMENA  AND  STABILITY  OF  PROBABILITY  DENSITIES 


I.  Introduction 

We  are  interested  in  the  overall  time  dependence  of  the  relaxation  of  a 

perturbed  entity  to  equilibrium.  The  overall  time  dependence  of  the  perturbed 

part  of  an  entity  is  frequently  tracked  in  terms  of  the  step  function  response 

function^  <j>(t).  The  function  <J> ( t )  has  maximal  value  at  t*0,  is  continuous,  and 

vanishes  as  t-*-«®.  The  function  4>(t)  nay  be  identified  as  the  positive  time  part 

of  an  even  function  $(t)  which  may  be  related^*^  to  the  dissipative  part  of  a 

scalar  impulse  response  function  x"(0  by 

X " ( t )  =  -  ,  —  <  t  <  «  (1.1) 

dt 

In  the  following  we  only  consider  situations  t*here 

lim  X”(t)  =■  0  (1.2) 

t*<» 

so  that  the  Fourier  transform  of  x'^t)*  namely  the  dissipative  part  of  the  sus- 
ceptibility  x"(w),  has  only  a  continuous  spectrum  with  no  singularity  at  u)»0. 
Further,  since  the  Fourier  transform  of  x"(b)  or  the  dissipative  part  of  the 
susceptibility  has  all  moments,  all  the  derivatives  of  x”(t)  and  hence  also 
$(t)  at  t*0  must  be  finite. 

The  dissipative  impulse  response  function  x"(t)  may  be  related  by  pertur¬ 
bation  theoretical  developments  such  as  those  of  Kubo,^  Case,^  and  Martin^  to 
the  equilibrium  average  of  a  two-time  scalar  commutator  of  a  configurational 
coordinate  operator  A(t), 

X-(t-t’)  -  i_  Tr{PofA(t),  A(t 1 )  ]  }  (1.3) 

ifi 

Here  p0  is  the  canonical  equilibrium  density  matrix  and  we  have  taken  A(t)  to 

be  the  perturbed  part  only.  Then  it  is  clear  that  x"(t)  and  also  $(t)  are 

functions  only  of  a  time  difference  t. 
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The  equation  for  $(t)  corresponding  to  x"^)  Eq.  Cl  -  3)  is2>^ 

*(t)  “  2  E  Wm~Wm*  |  <m  I A I  m '  >  | 2  cos(E0-Em.  )t/-Ti  (1.4) 

m>m'  Em.-Em 

Here  Em  is  an  eigenvalue  of  the  unperturbed  Hamiltonian,  Wm=exp(-BEm)/Z, 
Z”Eexp(-0Em)  and  B*l/kT.  It  is  clear  that  4>(t)  is  even  in  t.  In  the  high 
temperature  limit,  equivalent  to  a  classical  limit,  $(t)  is  /just  B/2  times  the 
canonical  average  of  the  anticommutator  (A(t),  A(0)}  which  may  be  identified 

with  a  real  autocorrelation  function.  In  any  event  $(t)  shares  with  the  auto¬ 
correlation  function  of  a  weakly  stationary  process  the  general  properties  that 
it  is  a  two-time  function  dependent  only  on  time  differences,  i.e.  location 
invariant,  and  continuous  for  all  time.  We  have  also  in  addition  that  $(t)  has 
other  temporal  behavior  determined  by  the  physical  properties  of  $(t)  mentioned 
above,  viz.  it  is  maximal  and  differentiable  to  all  orders  at  t=0,  and  mono- 
tonically  vanishing  as  t><».  Furthermore  from  the  causality  of  $(t), 
we  infer  x”(<*>)  is  one  of  a  Hilbert  transform  pair. 2  Consistency  with  this 
latter  property  of  x"^)  leads  us  to  the  requirement  that  $(t)  is  square 
integrable  for  -«<t<<»  i.e.  L2 (-«.«).  Also  from  Eq.  (1.4)  we  may  verify  that 
since  $(t)  is  continuous  at  t"0,  it  is  also  continuous  at  neighboring  points, 
e,  and  hence  for  all  t.  Moreover  if  X(t)  is  an  arbitrary  function  whose 
Fourier  transform  exists,  then  from  Eq.  (1.4)  it  follows  that 

//X(tH(t-t')X*(t,)dtdt'  =■  E  (Wm~Wm*  )  ]  <m  |  A  |  m  *  >  I  2  . 

m,tra’  Em.-En 

-  i  OVEm'H 

|/X(t)e  *  dt| 2  >  0  (1.5) 

3ochner's  theorem^  may  then  be  used  to  prove  that  <t(t)  is  proportional  to  a 
characteristic  function. 
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We  find  it  convenient  to  deal  with  the  normalized  function 
C(t|0)  =  *(t)/*(0)  (1.6) 

where  we  have  made  the  two  time  nature  of  C(t|0)  explicit.  By  its  definition 
C (0 | 0 )  -  1  so  that  it  is  in  fact  a  characteristic  function  and  we  shall  refer 
to  it  as  such  in  the  following. 

A  major  purpose  of  this  paper  is  to  construct  a  class  of  characteristic 
functions  with  behavior  that  is  consistent  with  the  mathematical  and  phenomeno¬ 
logical  constraints  on  the  step-response  function  <&(t).  These  requirements  are 
satisfied  by  allowing  a  continous  C(t|0)  to  be  only  differentiable  almost 
everywhere.  This  leads  to  a  piecewise  character  of  C ( t | 0 )  which  will  be  dis¬ 
cussed  in  detail  in  the  next  section.  The  resulting  construction  is  a  C(t|0) 
which  has  the  form  of  appropriate  exponentials  of  monomial  functions  of  t, 
namely  -( | t | termed  monomial  exponentials,  on  appropriate  non-overlapping 
segments  of  time.  For  example,  for  k*2  we  have  a  segment  of  quadratic  exponen¬ 
tial  behavior,  etc.  Continuity  conditions  at  the  joining  points  of  the  seg¬ 
ments  provide  relations  among  the  t^.  It  is  pointed  out  that  a  C(t|0)  so 
defined  is  in  fact  the  Fourier  transform  of  a  probability  density. 

With  a  probability  density  in  hand,  we  turn  to  another  major  purpose  of 
the  paper  in  section  III  and  consider  the  description  of  a  total  system  built 
up  from  subsystems  that  share  the  same  generic  temporal  behavior.  We  find  the 
conditions  under  which  the  probability  density  for  the  total  system  also  has 
the  same  piecewise  characteristic  function  as  that  of  the  subsystems.  A  new 
type  of  stability  for  probability  densities  is  thereby  obtained. 

In  section  IV,  we  compare  the  general  properties  of  the  piecewise  C ( t | 0 ) 
for  t>0  with  the  step-response  function  in  a  model  of  relaxation  in  complex 
systems  previously  introduced^  and  discussed^  by  us.  The  latter  in  fact  shares 
two  important  features  with  the  C(t|0)  constructed  in  section  III.  Namely  it 
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has:  (1)  a  long  time  fractional  exponential  regime  preceded  by  a  linear  expo¬ 
nential  regime;  and  (2)  a  relationship  between  the  relaxation  times  or 
r-parameters  in  the  two  regimes.  The  model  is  reviewed  briefly  as  a  way  to 
gain  insight  into  the  meaning  of  the  piecewise  C(t|0)  and  the  new  type  of  sta¬ 
bility  for  probability  densities  discussed  in  section  IV. 

II.  Temporal  Behavior  of  the  Characteristic  Function 

We  want  now  to  consider  the  consequences  of  the  various  phenomenological 
and  mathematical  constraints  on  the  temporal  behavior  of  the  characteristic 
function  C(t|0)  introduced  in  section  I.  The  properties  already  indicated  in 
that  section  are  the  general  ones  of  location  invariance  and  continuity,  and 
the  more  special  ones  of  evenness  in  t,  differentiability  to  all  orders  at  t*0 
(from  the  existence  of  all  moments  of  the  dissipative  part  of  the  suscepti¬ 
bility),  maximum  value  at  t=0,  and 

It  should  be  noted  that  the  physics  of  all  our  considerations  is  nonrela- 
tivistic  so  the  details  of  behavior  at  or  near  t=0  is  not  actually  included  in 
our  formulation.  The  best  we  can  hope  for  near  t*0  is  a  self-consistent  phe¬ 
nomenology.  However,  this  self-consistency  leads  to  some  interesting  con¬ 
straints.  From  the  differentiability  property  of  C(tlO)  at  t*0,  linear  (and 
all  other  odd)  terms  in  its  power  series  expansion  around  t=*0  must  vanish.  In 
particular,  C(tlO)  cannot  then  be  taken  to  have  a  linear  exponential  form  at 
ta0.  An  even  more  stringent  condition  arises  from  the  fact  that  all  the  deriv¬ 
atives  of  C(t|0)  at  t=0  must  be  finite.  Further  as  already  noted,  from  the 
physical  meaning  of  step  response,  C(t|0)  is  also  maximal  at  t”0.  A  self- 
consistent  (but  not  unique)  choice  for  C ( t | 0 )  at  or  near  t*0  is  a  quadratic 
exponential  or  a  Gaussian  form  centered  on  t=0. 


For  values  of  tX) ,  we  expect  from  empirical  evidence  covering  most  parti¬ 
cular  cases  that  C(t|0)  mav  appear  to  decay  in  one  of  three  general  forms:  (a) 
quadratic  exponential  (Gaussian);  (b)  linear  exponential  (most  usually  calcu¬ 
lated,  and  most  often  observed  in  simple  systems);  (c)  fractional  exponential 
(for  complex  systems).  Further  for  very  long  times,  the  real  L^(-«,<»)  charac¬ 
teristic  function  considered  here  is  constrained  to  be  consistent  with  the 
Paley-Wiener  bound® 

Mm  C(t|0)  >  exp-(  1 1 1 /Ta)  ,  0<ct<l  (2.1) 

1 1 1  -*■« 

where  Ta  is  a  positive  real  scale  parameter.  It  may  be  noted  that  the  Paley- 
Wiener  bound  is  a  fractional  exponential  so  the  latter  is  the  fastest  allowable 
decay  as  1 1 1  -*•<». 

Thus  for  any  C(t|0)  there  are  conditions  near  or  at  t=0  and  for  very  long 
times  that  must  be  satisfied.  At  intermediate  times,  there  may  be  physically 
determined  behavior  that  differs  from  at  least  one^  and  possibly  both  of  the 
behaviors  respectively  around  t=0  and  for  very  long  times.  We  now  consider  a 
generic  C ( t | 0 )  that  incoroorates  all  these  features  and,  dependent  upon  the 
choice  of  parameterization,  can  be  used  to  describe  any  of  the  observed  mono¬ 
mial  exponential  time  behaviors. 

We  consider  C(t|0)  to  be  piecewise  continuous  with  appropriate  monomial 
exponential  behavior  in  given  nonoverlapping  segments  which  cover  -®<t<»  com¬ 
pletely.  The  conditions  around  t=0  and  for  a  physical  decay  regime  that  has  a 
quadratic  monomial  behavior  can  be  satisfied  by  taking  C(t[0)  to  have  the  form 
exp-(t/i2)^  in  an  initial  segment  pair  .-joined  at  t=*0,  i.e.  for  0<  1 1 1  <Tq.  Here 
t2  is  a  positive  real  scale  parameter  and  Tq  is  also  nositive  real.  The  pos¬ 
sibility  of  a  linear  exponential  regime  that  often  occurs  empirically  and  is 
the  form  obtained  in  standard  decay  calculations  can  be  allowed  by  choosing 
C(t|0)  to  have  the  form  exp-ltl/ti  in  a  second  segment  pair  Tq<  f  1 1 <TL.  Here 
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gain  n  is  a  positive  real  scale  parameter  and  “l  is  positive  real.  To  allow 
or  complex  physical  systems  that  have  a  fractional  exponential  decay  and  for 
ong  time  behavior  at_  the  Paley-Wiener  bound,  C(t|0)  is  taken  to  have  the  form 
xp-( I t | /xa)a,  0<a<l  for  Tj_i<  1 1 1  <®.  Again  Tl  is  positive  real.  The  continuity 
f  C(t|0)  is  assured  by  continuity  conditions  at  its  -joining  points  ±Tq,  ±Tl. 
he  time  behavior  of  our  generic  C(t|0)  may  be  summarized  as  follows. 

f  exp-(  t/t2)2  ,  -T0<t<TQ 

exp-(Tq/x2)2  =  exp-( Tq/x^ )  ,  It!  =  Tq 

(t|0,  T0,  Tl;  xk)  =»  J  exp-|t|/Ti  ,  Tq< | t ! <TL 

\  exp-(TL/xi)=exp-(TL/Ta)a  ,  0<a<l ,  1 1 1  =  TL 


exp-( I t | /xa)a,  0<a<l,  T^<  1 1 1  <® 


(2.2) 


The  symbol  xk  in  C(t|0,  Tq,  T^;  xk)  represents  any  of  the  k=a,  1,  2  since 
nly  one  of  the  xk  is  an  independent  parameter.  From  the  continuity  condi - 
ions,  it  follows  that 


2  *  V1 


(2.3a) 


ta  »  Tl01-1!! 
a 


(2.3b) 


he  choice  of  xk  as  a  label  will  be  made  on  physical  grounds  as  described  in 
ore  detail  below.  For  the  moment  we  write 


C(t|0,  T0,  Tl;  tk)  =  CT  (tlO) 

k 


(2.4) 


nd  choose  tk  to  be  that  xk  which  falls  within  the  positive  t  segment  with  k 


xponential  behavior.  This  can  happen  for  only  one  of  the  k  values.  For 

xample,  if  Tq<x^<T^  we  would  write  CT^(t|0)  and  note  that  x2  =  Tqxi>T2, 

a  =.  T«~1t  ^<Ta  so  that  x2  is  not  in  [0,Tq]  and  xa  is  not  in  [Tl,®].  Similarly 
a  L  L 

or  CT  (t|0),  x2  is  in  [0,Tq]  but  xj  is  not  in  [Tq,Tl1  and  xa  is  not  in  [Tl,®]. 


Finally  for  CT  (t|0),  ta  is  in  [TL,®)  but  T2  is  not  in  [0,Tq]  and  is  not 
a 

in  [Tq.TjJ.  We  shall  see  that  CT^(t|0)  is,  for  centrally  placed  on  a 
logarithmic  time  scale,  dominated  by  its  behavior  in  the  region  in  which 
falls.  We  shall  call  CT^(tjO)  a  k-dorainant  characteristic  function  in  such 
cases.  However,  before  addressing  this  particular  point,  we  want  to  note  some 
general  properties  of  CT^(t|0). 

First  it  is  clear  that  although  CT  (t|0)  is  continuous,  it  is  only  piece- 

k 

wise  differentiable  with  discontinuities  in  derivative  at  t=±TQ,  iT^.  Further, 
location  invariance  is  maintained  by  the  requirement  that  a  translation  of  the 
origin,  e.g.  by  6,  is  accompanied  by  a  like  translation  of  time  segment  end¬ 
points  e.g.  -T^>- Tl+6  and  Tl+Tl+S,  etc.  Also  we  note  that  CT  (tlO)  is  the 

k 

Fourier  transform  of  a  probability  density.  This  follows  because  the  nature  of 
the  probability  distribution  function  associated  with  the  characteristic  func¬ 
tion  depends  on  the  manner  in  which  the  latter  behaves  as  1 1 1 If  it  goes 
to  zero  as  with  our  (-«,<■)  class  of  functions,  then  a  probability  distribu¬ 
tion  function  is  continuous.  For  completeness  it  may  be  noted  that  if  the 
characteristic  function  tends  to  a  constant,  D,  with  0<D<1 ,  the  associated 
probability  distribution  function  is  singular  and  may  be  of  the  form  of  a 
Cantor  set,  for  example.  If  it  tends  to  unity,  then  the  probability  distribu- 
i on  function  is  discrete. 

We  return  new  to  the  discussion  of  the  generic  character  of  Cx  (t|0) 

k 

which  is  readily  exhibited  by  considering  three  mutually  exclusive  cases  of 
parameter  magnitude:  (a)  0<<T2<<T^;  (b)  Tq<<X[<<T{j5  (c)  TL<<xa<<®,  Case  (a) 

corresponds  to  2-dominant  behavior.  The  malor  part  of  the  decav  of  C _  (tlO) 

L  *1 

appears  as  a  quadratic  exponential,  and  the  linear  and  fractional  exponential 
parts  are  submerged  in  a  long  time  tail  of  small  magnitude.  On  the  basis  of 
Its  direct  effect,  we  could  here  omit  the  linear  exponential  part.  However, 


present  construction  requires  the  continuity  condition  Eq.  (2.3a)  that  pro- 
es  a  relation  between  12  and  tj.  Thus  the  linear  exponential  region  would 
revealed  by  comparing  functional  relations  for  a  measured  T2  with  corre- 
nding  functional  relations  for  a  calculated  tj.  Finally,  the  fractional 
onential  or  slower  behavior  must  be  included  at  long  times  to  meet  the  con- 
ion  of  the  Palev-Wiener  bound. 

Case  (b)  corresponds  to  1-dominant  behavior.  The  major  part  of  the  decay 
CT  ( t | 0 )  appears  as  a  linear  exponential.  The  latter  is  of  course  the  most 
monly  observed  behavior  in  simple  systems  and  is  the  usual  result  of  theo- 
ical  calculations  based  on  a  random  phase  approximation  or  the  introduction 
a  phenomenological  friction  terra.  The  quadratic  exponential  segment  of 
(t|0)  is  submerged  in  small  deviations  at  times  short  compared  to  T}  but 
or  some  equivalent  must  be  included  to  meet  the  conditions  on  CT  (t|0) 
t=0.  Also  again  a  fractional  exponential  or  slower  time  behavior  must  be 
luded  to  meet  the  Paley-Wiener  bound  at  long  times  even  though  that  region 
a  small  magnitude. 

Finally,  case  (c)  corresponds  to  a-dominant  behavior.  The  major  part  of 

decay  of  CT  (t|0)  appears  as  a  fractional  exponential.  The  quadratic  expo- 
a 

tial  segment  has  an  immeasurable  effect  at  short  times  but  it  or  its  eauiva- 
t  is  essential  to  meet  conditions  at  t=0.  The  linear  exponential  regime 
e  provides  only  small  deviations  that  would  most  probablv  also  not  be  noted 
measurements.  Thus  again  as  in  case  (a)  the  linear  exponential  segment 
ht  at  first  sight  appear  to  be  unnecessary.  However,  if  there  is  a  linear 
onential  regime,  the  measured  parameter  xa  can  be  compared  to  a  calculated 
bv  means  of  the  continuity  condition  Eq.  (2.3b).  This  continuity  condition 
n  then  provides  an  important  relation  that  is  subject  to  verification.  In 
t,  a  similar  relation  between  relaxation  time  parameters  respectively  in  the 


linear  and  fractional  exponential  regimes  is  a  key  feature  of  the  model  of 
relaxation  in  complex  systems^#?  to  be  discussed  below  in  section  IV. 

We  can  gain  further  insight  into  the  time  behavior  of  CT^(t|0)  by  con¬ 
sideration  of  the  instantaneous  decay  rate  W_  (t|0)  which  is  defined  to  be 

k 


WT  (t | 0)  =  -(d  CT  (t|0)/dt)/CT  (t|0) 
k  k  k 


=  -  d  £nCT  (t|0)/dt 
Tk 


t  *  ±  tq>  Tl 


(2.4) 


We  designate  WT  (t|0)  to  be  W2(t|0),  W^(t|0),  and  W  ( 1 1  0 )  respectively  in  the 
k 

segments  where  CT  (t|0)  has  quadratic,  linear  and  fractional  exponential  behav- 
k 

vior.  Since  CT  (t|0)  is  even,  it  is  sufficient  to  consider  t>0.  In  the  quad- 
k 

ratic  exponential  segment.  Wo ( t I 0)»2 t/t^  which  approaches  a  maximum  value 

2 

W2(To|0)»2Tq/t^**2/t[  as  t  approaches  Tq.  On  the  other  hand  Wj  ( 1 1  0)*1  /x\  for 

all  Tq<t<TL  so  there  is  a  discontinuity  in  WT^  (t|0)  at  t=Tq  where  it  falls  to 

one  half  its  limiting  maximum  value  in  the  quadratic  segment.  Similarly  there 

is  a  discontinuity  at  t*Tq  for  Wa(t|0)  *  ata“^/xa  has  a  maximum  limiting  value 

a 

Wa(TLIO)  ■  a/xy  <  1  / t i .  Ir  is  clear  that  this  sequential  behavior  for  W_  (t|0) 

k 

is  independent  of  the  choice  of  magnitude  t's  and  so  occurs  for  all  possibili¬ 
ties  of  the  k  in  k-domi nance. 

The  sharpness  of  the  discontinuities  in  WT  (t|0)  is  an  artifact  of  our 

piecewise  construction  for  CT  (t|0).  Nevertheless  decay  rates  depend  on  system 

k 

interactions  so  these  or  perhaps  somewhat  more  smooth  behavior  can  be  given  a 
physical  meaning.  They  can  be  viewed  as  due  to  changes  in  the  nature  of  the 
interactions  in  the  system  at  |t!=Tq,  Tl*  Such  changes  can  be  thought  of  as 
temporal  transformations  of  the  environment  in  which  the  relaxation  is 
occurring.  These  transformations  can  be  viewed  as  the  analogs  of  thermodynamic 


phase  transformations. 
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The  CT  (t 1 0)  of  Eq.  (2.2)  could  be  readily  generalized  by  addition  of 
Tk 

other  monomial  exponential  intermediate  segments  or  appropriate  short  time  or 

long  time  segments.  Such  a  generalization  would  also  involve  corresponding 

additions  (and  discontinuities)  in  WT  ( t | 0) .  The  point  we  wish  to  emphasize 

k 

however  is  that  the  generic  CT^(t|0)  defined  in  Eq.  (2.2)  satisfies  all  the 

mathematical  requirements  for,  and  can  provide  an  adequate  description  of  the 

most  usually  phenomenological  behaviors  of  real  step-response  functions.  A 

simple  and  useful  property  of  CT  (t|0)  so  constructed  is  that  it  has  monomial 

k 

exponential  behavior  in  all  segments.  This  has  provided  simple  continuity  con¬ 
ditions,  Eqs.  (2.3a,b),  and  further,  provides  the  basis  for  consideration  in 

the  next  section  of  stability  properties  for  the  probability  density  PT  (id) 

k 

belonging  to  CT  (t|0). 

k 


III.  k-dominant  Stability  for  Probability  Densities 

We  turn  now  to  a  discussion  of  the  stability  properties  of  PT^(a>).  The 

probability  density  PT  (ui)  corresponding  to  CT  (t|0)  is  a  function  of  id  as  a 

k  k 

random  variable  that  incorporates  properties  that  reflect  the  particular 


k-dominant  and  piecewise  properties  of  CT  (t|0). 

k 

CT  (t|0)  =  /coPT  ( ui)  cosut  dco 
k  ~~  k 

P x  (u» )  =  /*  CT  ( 1 1  0 )  cos  u>t  dt 

k  2-rr  —  k 


(3.1a) 

(3.1b) 


We  note  that  CT  ( 1 1 0 )  remains  invariant  under  scaling  of  tir..e  bv  a  numerical 
Tk 

factor.  From  Eqs.  (2.2)  and  (2. A), 

(3.2) 

This  property  merely  reflects  the  fact  that  the  choice  of  scale  in  temporal 


C(At | 0 ,  ATq,  ATl;  Atk)  =  CAt  (At|0)  =  CT  (t|0) 

k  k 


measurements  does  not  affect  the  behavior  of  CT  (t|0).  On  the  other  hand  from 

Tk 

Eq.  (3.1a),  the  probability  density  must  then  obey  the  following  scaling  rela¬ 
tion 
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(3.3) 


pAt  ”  APt  (“A) 
k  k 


Consider  now  N  subsystems  each  of  which  is  characterized  by  a  member  of 
the  set  of  independent  and  identically  distributed  stochastic  variables  We 

want  then  a  total  system  made  up  of  these  subsystems  and  characterized  by  a 
normalized  summand  variable 


-1  N 

a)  =*  B  £  u)j  (3.4) 

Nk  i=1 


such  that  the  probability  density  for  oj  is  also  the  Fourier  transform  of  the 
same  k-doninant  characteristic  function  as  that  for  the  individual  systems. 
The  normalization  function  is  expected  to  depend  on  both  N  and  k.  Such  a 
probability  density  will  be  called  k-doninant  stable.  We  write 

_1  N  N 

PT  (u>)  -  <5(B  £  iiii-(o)  H  PT  (o»1)dwi  ,  (3.5) 

k  -*  -*  Nk  k 


which  from  Eq.  (3.3)  can  also  be  written  in  the  form 
N  N  _  N  _ 

PT  (u>)  «  B  6(  £  u  -w)  n  PT  (B^  (Di  )du)i  (3.6) 

k  l-«  ■«  mi  i=i  k 


Here 

“i 


-1 

B  uj-jl 
Nk 


(3.7) 


Transforming  the  <5-function  to  t-space  and  carrying  out  the  N  uif  integrals, 
we  obtain 


P f  (ui)  -  /-  tCB  T  (t|  0)  ]N  cos  u)t  dt 

k  2tt  “c#  r,k  k 


(3.8) 


The  time  arguments  in  the  products  of  Cg  T  (t|0)  in  the  integrand  are  the  same 

Nk  k 

so  the  resulting  product  of  N  Cg  ,  (t|0)’s  has  the  piecewise  structure  of 

Nk  k 

Co  T  (t|0).  However,  Eq.  (3.8)  is  not  a  proper  equation  unless 
Nk  k 

[Cg  T  <t 1 0) ]N  -  cT  (t| 0)  (3.9) 

Nk  k  k 

Using  the  monomial  exponential  property  of  all  the  segments  of  Cg  T  (t|0)  and 

Nk  k 

continuity  conditions  Eqs.  (2.3a,b)  as  appropriate,  we  find  Eq.  (3.9)  follows 
if 

%k  =  Nl/k  (3.10) 

We  have  the  result  then  that  k-dominant  stability  is  achieved  for  a  prob¬ 
ability  density  of  a  summand  variable  that  is  normalized  by  N* /k.  For  k*2 , 
N^/2  provides  an  enhancement  in  ui  with  respect  to  the  id^.  For  k*l ,  tu  and  the 
individual  oi^  coincide  while  for  k=a,  0<a<l ,  oj  is  reduced  with  respect  to  an 
individual  to^. 

We  can  view  the  normalization  as  a  manifestation  of  an  interaction  envi¬ 
ronment  that  is  introduced  when  the  independent  subsystems  are  put  together  to 
form  the  total  system.  However,  once  the  total  system  is  formed,  it  can  be 
viewed  as  a  building  block  subsystem  for  a  larger  total  system.  It  follows 
that  the  interaction  environment  for  the  single  system  must  be  self-consistent 
with  the  interaction  environment  for  the  total  system.  An  important  question 
is  how  the  interaction  environment  is  first  set  up  on  a  microscopic  basis  or 
alternatively,  if  k-dorainant  stability  is  taken  as  a  primary  principle,  how  a 
microscopic  interaction  environment  is  arranged  to  form  the  appropriate  build¬ 
ing  block.  In  the  next  section  we  consider  a  model  of  relaxation  in  complex 
systems  that  can  be  viewed  as  providing  a  microscopic  basis  for  a-dorainant 
behavior.  If  k-dominant  stability  is  taken  as  the  primary  principle,  all 
models  of  relaxation  must  be  constructed  in  such  a  way  to  yield  k-dominant 


behavior 


In  general,  this  complementarity  of  microscopic  basis  and  k-dominant 


stability  may  be  viewed  in  the  following  physical  way.  Theoretical  descriptions 
of  relaxations  in  a  complex  physical  systems  such  as  glasses  and  polymers  are 
usually  based  on  a  microscopic  relaxation  species.  The  latter  can  e.g.  be  a 
dipolar  group,  a  hopping  charge,  a  spin,  or  an  entire  polymer  chain  or  segment 
of  it.  A  particular  relaxation  species  i  in  a  complex  environment  and  inter¬ 
acting  with  a  nominal  heat  bath,  will  have  behavior  described  in  terms  of  a 

.  i  . 

a  relaxation  function  ^1(t),  or  a  characteristic  function  CT  ( 1 1  0 )  or  a  proba- 

i  k 

ability  density  PT  (oi).  Measurements  in  relaxation  are  macroscopic  quantities 
k 

such  as  transient  current,  stress,  magnetization,  volume  and  shear  flow.  The 
* 

stochastic  variable  oi  of  the  macroscopic  quantity  is  a  normalized  sum  of  the 

form  of  Eq.  (3.4).  It  is  generally  assumed  in  all  considerations  of  relaxation 

that  the  macroscopic  <j>(t)  is  identical  to  the  microscopic  <)>*•(  t).  This  assump- 

i 


tion  is  equivalent  to  the  assumption  that  PT^(u)  is  a  stable  probability  den¬ 


sity  or  in  our  description,  k-dominant  stable.  From  the  results  of  this  sect¬ 


ion,  these  equivalent  assumptions  are  consistent  with  the  requirement  that 


$*(t)  be  the  positive  time  part  of  a  k-dominant  CT  ( 1 1 0)  of  the  form  given  in 

k 


Eq.  (2.2). 

As  a  final  comment  in  this  section,  it  is  interesting  to  consider  the 
relationship  between  Levy  stability^*  ^  and  the  k-dorainant  stability  intro¬ 
duced  here.  While  we  have  restricted  ourselves  to  real  characteristic  func¬ 
tions  that  must  be  piecewise  with  multiple  monomial  exponential  behaviors. 


Levy  deals  with  characteristic  functions  with  magnitude  equal  to  a  single  mono¬ 
mial  exponential  with  monomial  power  g,  0<g<2.  Except  for  the  possibility  of 
the  quadratic  exponential  behavior,  such  monolithic  behavior  cannot  meet  the 
physical  conditions  on  the  autocorrelation  function  at  t*0.  For  a  real  charac¬ 
teristic  function,  the  monolithic  quadratic  exponential  is  ruled  out  by  the 


physical  unrealistic  requirement  that  its  instantaneous  decay  rate  is  propor¬ 
tional  to  t  and  therefore  blcws  up  as  t-*-00  and  such  behavior  is  never  observed. 
For  any  characteristic  function,  the  limitation  to  a  quadratic  exponential  mag¬ 
nitude  at  intermediate  times  is  not  sufficiently  general  to  provide  a  satisfac¬ 
tory  basis  for  the  description  of  decay  phenomena.  Thus  from  the  viewpoint  of 
physics,  the  present  generalization  to  k-dominant  stability  and  its  concomitant 
piecewise  characteristic  function,  or  some  further  generalization  yet  to  be 
made,  is  a  necessary  step  beyond  Levy  stability. 

IV.  Comparison  with  the  Ngai  Model  of  Relaxation 

We  now  take  advantage  of  the  fact  that  CT^(t|0)  for  t>0  can  be  identified 
as  a  step  response  function  to  compare  our  constructed  CT^(t|0)  with  the  relax¬ 
ation  function  obtained  in  the  so-called  Ngai  model  of  relaxation.  The  latter 

model  predicts:  (1)  a  fractional  exponential  decay  in  a  time  regime  t>>o)“^ 

c 

when  u)cTa»l ;  (2)  a  relationship  between  a  relaxation  time  in  a  linear  expo¬ 
nential  regime  t<<oj~1  and  xa.  Here  u)c  is  a  cutoff  frequency  that  is  a  param- 

c 

eter  of  the  model  as  described  below  and  we  have  taken  the  liberty  of  using 

tu,  k=*a,I  to  suggest  the  direct  connection  with  our  CT  (t|0).  It  is  clear 

k 

that  the  model  predictions  are  consistent  with  k-dominant  behavior  for 

CT  (t|0)  and,  as  we  shall  see  in  detail  below,  the  continuity  condition  Eq. 
k 

(2.3b).  Thus  the  model  may  indeed  provide  insight  into  the  physics  of  the  for¬ 
mation  of  CT  ( 1 1  0 )  and  the  question  of  the  microscopic  origin  of  the  building 
k 

blocks  of  k-dominant  stability. 

To  provide  context  for  such  a  discussion,  we  briefly  review  the  concepts 
in  the  Ngai  model.  The  model  considers  a  relaxing  entity  that  is  coupled  to  a 
heat  bath  on  a  microscopic  level.  In  accordance  with  the  conventional  descrip¬ 
tion  of  relaxation,  after  some  induction  time,  say  t>Tq,  the  system  initially 


,  -  .  * 
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relaxes  according  to  a  linear  exponential.  The  complexity  of  the  system  is 
modelled  in  terms  of  correlations  imposed  on  the  simple  system.  These  correla¬ 
tions  are  represented  by  modes  which  are  distinct  from  the  heat  bath  but  which 
also  interact  with  the  relaxing  entity.  These  modes  were  called  correlated 
state  excitations  in  Ref.  6.  The  correlated  state  excitations  themselves 
modell,  in  terms  of  the  energy  level  spacings  of  the  Gaussian  Orthogonal 
Ensemble  (GOE)  of  random  matrix  Hamiltonian  theory.^  The  level  spacing  dis¬ 
tribution  of  the  GOE  gives  the  energy  density  function  N(E)  of  the  correlated 
state  excitations  which  are  also  called  level  spacing  excitations. ^  N(E)  varies 
linearly  with  E  up  to  an  upper  cutoff  Ec*tlajc  that  defines  the  low  frequency 
regime,  levels  off  to  a  plateau  at  intermediate  energies,  and  vanishes  (linear) 
exponentially  at  high  energies.  As  noted  above,  long  and  short  time  regimes 
for  the  model  are  quantified  in  terms  of  mc,  i.e.  ioct>>l  means  t  in  the  long 
time  regime  while  u>ct<<l  means  t  is  in  the  short  time  regime  (for  us,  inter¬ 
mediate  time  regime). 

The  relaxing  entity  is  then  not  only  coupled  to  the  heat  bath  but  also  to 
the  energy  level  spacing  or  correlated  state  excitations.  The  latter  coupling 
strength  is  assumed  to  be  independent  of  the  energy  of  the  excitations.  It  is 
important  to  emphasize  that  the  level  spacing  excitations  are  not  part  of  the 
heat  bath.  The  excitations  are  in  equilibrium  prior  to  system  excitation,  are 
driven  out  of  equilibrium  during  the  relaxation  by  their  interaction  with  the 
measured  relaxing  entity,  and  eventually  return  to  equilibrium  in  a  process  not 
addressed  specifically  in  the  model. 

The  interaction  between  the  relaxing  entity  and  the  level  spacing  excita¬ 
tions  leads  to  a  long  time  fractional  exponential  decay  as  a  manifestation  of 
the  modification  of  the  instantaneous  decay  rate  from  W^(t)“T~^  for  u>ct<<! 
Wa(t)  -  Wj(t)  exp[(a-l)y]  (a>ct)a-1  ,  0<a<l ,  u>ct»l  .  (4.1) 


Here  y  is  the  Euler  constant.  Instead  of  a,  the  model  is  usually  formulated  in 
terms  of  the  parameter 


n  *  1-a  (4.2) 

where  n  *  2C|v|2  .  (4.3) 

Here  C  is  the  slope  of  the  linear  portion  of  the  density  of  states  N(E)  and 
| V | 2  is  the  averaged  squared  magnitude  of  the  coupling  interaction  between  the 
relaxing  entity  and  the  level  spacing  excitations.  As  mentioned  above,  |V|  is 
assumed  to  be  independent  of  the  energy  of  the  excitation. 

The  long  time  step  response  function  is  written  as 

4>(t)  -  exp-(t/Ta)1-n  ,  a)ct»l  (4.4) 

It  should  also  be  noted  that  Eq.  (4.4)  applies  for  experimental  observations 

where  oJcra»l  so  <p ( t )  meets  the  condition  (c)  of  ci-dominance.  Given  Eq.  (4.4), 

the  instantaneous  decay  rate  for  aict»l  is 

Wa(t)  *  ata“^ / Ta  (4.5) 

a 

Consistency  of  Eqs.  (4.1)  and  (4.5)  requires  that 

xa  *  a  exp[(l-a)y]  (4.6) 

a  c 

Equations  (4.4)  and  (4.6)  are  the  predictions  of  the  Ngai  model  which  are 

readily  susceptible  to  experimental  verification.  For  example,  let  us  consider 

the  effect  of  elevated  temperature  on  relaxations  in  complex  systems  such  as 

amorphous  polymers  and  glasses  which  at  normal  temperatures  have  the  a-dominant 

form  of  Eq.  (4.4)  and  xa  is  thermally  activated.  The  Ngai  model  has  an  addl- 

^ional  prediction  given  by  Eq.  (4.6)  which  relates  xa  to  tj.  If  xa  has  an 

Arrhenius  temperature  dependence  with  activation  energy  E*,  then  xj  will  have  a 

a 

different  but  predictable  activation  energy  of  Ea  =  aE*.  More  generally,  a 

a 

dependence  f($)  of  ra  on  a  physical  parameter  5,  such  as  molecular  weight  M  of 


-I  Z? 


a  polymer  chain  in  polymer  melt  relaxation  or  isotope  mass  m  of  alkali  ion  in 
conductivity  relaxation  of  glasses,  will  correspond  to  a  different  lf(c)]a 
dependence  on  £  for  x j .  These  predictions  have  been  verified  in  many  different 
complex  systems.^ 

An  interesting  possibility  that  is  allowed  in  the  Ngai  model  is  the 

cross-over  from  a-dominant  to  1-dominant  behavior  by  changing  the  conditions  of 

physical  measurement.  To  discuss  this  point  further,  we  complete  the  contact 

between  the  Ngai  model  and  the  construction  of  CT  (t|0)  by  identifying  T^  in 

a 

Eq.  (2.3b)  to  be 

Tl  *  al/(a“l)exp (— y)u)“l  (4.7) 

c 

Then  we  can  consider  a  change  in  physical  conditions  by  which  Tl  does  not  vary 
significantly  but  where  ta  is  changed  by  many  orders  of  magnitude  so  that  a 
crossover  occurs.  In  such  a  case,  the  new  x^  can  obey 

?1  «  TL  (4.8) 

and  the  behavior  is  1-dominarit  because  it  is  now  rj  that  is  in  the  interval 
(Tq,Tl).  Perhaps  surprisingly,  such  a  situation  occurs  in  a  straightforward 
way  in  the  temperature  dependence  of  the  relaxation  of  several  materials 
including  the  glass-forming  fused  salt^  0. 4Ca(NC>3 )2*  (0. 6)KN03 ,  the  ionic  con¬ 
ductor  sodium  0-alumina,  ^  and  in  polymeric  systems.^ 

It  may  be  noted  that  most  presently  proposed  models^  of  relaxation 
encompass  neither  the  prediction  of  Eq.  (4.6)  nor  the  possibility  of  crossover 
of  type  of  dominant  behavior.  Although  the  Ngai  model  does  not  require  the 
sudden  discontinuity  by  changing  the  decay  rate  from  Wi=x“*  to  Wa(t)  at  the 
single  point  t*TL,  it  is  compatible  with  such  an  interpretation.  Then  instead 
of  the  more  abstract  discussion  of  the  temporal  analog  of  thermodynamic  phase 
transitions  in  section  II,  we  may  now  consider  Wj(t)  to  be  suddenly  dressed  by 
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means  of  interaction  with  correlation  or  level  spacing  excitations  at  to 
yield  the  decay  rate  Wa(t). 

Correlation  excitations  may  also  be  viewed  as  providing  the  vehicle  by 
which  a-dominant  stability  is  introduced  microscopically.  Thus  the  Ngai  model 
may  be  viewed  as  providing  the  first  building  block.  On  the  other  hand,  as 
indicated  in  section  III,  we  may  view  a-dominant  stability  as  a  primary  prin¬ 
ciple.  Then  the  microscopic  interactions  that  enter  the  first  building  block 
must  be  compatible  with  a-dominant  (or  more  generally  any  k-dominant)  stabil¬ 
ity. 
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